Geometric quantum gates robust against stochastic control errors 
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We analyze a scheme for quantum computation where quantum gates can be continuously changed 
from standard dynamic gates to purely geometric ones. These gates are enacted by controlling a set 
of parameters that are subject to unwanted stochastic fluctuations. This kind of noise results in a 
departure from the ideal case that can be quantified by a gate fidelity. We find that the maximum 
of this fidelity corresponds to quantum gates with a vanishing dynamical phase. 
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An essential prerequisite for quantum computation 
is the ability of maintaining quantum coherence and 
quantum entanglement in a information-processing sys- 
tem Unfortunately, since both these properties are 
very fragile against control errors as well as against 
unwanted couplings with environment, this goal is ex- 
tremely hard to achieve. To this end several strate- 
gies have been developed, most notably: quantum er- 
ror correction|2j,decoherence-free subspace|2(, and bang- 
bang techniquesjj]. 

Quantum computation implemented by geometric 
phases [EE10 is believed to be another approach which 
can be used to overcome certain kinds of errors H, 0, 0, 
ITU IT^ . fl3L IT3. Hlj . It has been shown a universal set 
of quantum gates can be realized by geometric phases. 
However, the statement that quantum gates achieved by 
this way may have built-in fault-tolerant features (due 
to the fact that geometric phases depend only on some 
global geometric properties) has still the status of a con- 
jecture and it has been the subject of some debate in the 
literature. Indeed this alleged resilience against errors of 
geometrical gates has been doubted by some numerical 
calculations with certain decohering mechanisms flfi Il7| . 
On the other hand, analytical results show that adiabatic 
Berry's phase itself may be robust against dephasin g|l8l 
and stochastic fluctuations of control parameters [19j. 
These latter provide a sort of indirect evidence of the ro- 
bustness of adiabatic geometric quantum computation; 
a more direct evidence would be given by a comparison 
between the fidelity of geometric gates and standard dy- 
namic gates in presence of an error source. So far, this 
kind of convincing evidence, clearly showing that geomet- 
ric quantum computation is robust against some realistic 
noise sources, is still missing. 

In this paper we shall consider a parametric family of 
quantum gates subject to stochastic fluctuations of the 
control parameters. The departure from the ideal (i.e., 
no-fluctuation) case can be quantified by gate fidelity. We 
compute such a fidelity in a model where quantum gates 
can continuously change from standard dynamic gates to 
purely geometric ones. We find that the maximum of 
fidelity corresponds to those cases in which the dynam- 
ical phase accumulated over the gate operation is zero. 



This provides a clear evidence of the robustness of nona- 
diabatic geometric computation in this specific scheme. 
The predictions presented here may be experimentally 
tested in some quantum computer prototypes. Moreover, 
this robustness of nonadiabatic geometric computation 
is relevant to experiments on quantum information pro- 
cessing; indeed it is generally believed that the control 
parameter nonuniformity is one of the most dangerous 
sources of errors for qubits in solid-state systems [20| . 
NMR H3 and trapped ions H3, etc.. 

Quantum computation via a pair of orthogonal cyclic 
states- For universal quantum computation, it is suffi- 
cient to enact two of noncommuting single-qubit gates 
and one nontrivial two-qubit gate. Before to study the 
fidelity of a quantum gate subject to noise, we recall a 
scheme to implement a universal set of quantum gates by 
using a pair of orthogonal cyclic states[I2, 03 ■ A single 
qubit gate given by 

_ f e 11 cos 2 ^ + e~ 11 sin 2 ^ i sin % sin 7 

\ i sin x sin 7 e 17 sin 2 ^ + e~ n cos 2 ^ 

can be obtained when a pair of cyclic states |±) can 
be found for a unitary operation i.e., E/W|±) — 

e ±l7 |±). Here x is related to the initial cyclic states, and 
7 is a phase accumulated in the gate evolution. Usually 
the total phase 7 consists of both geometric (7^) and dy- 
namic components (7^), and U is specified as a geometric 
gate if 7 is a pure geometric phase. Moreover, a condi- 
tional two-qubit gate can also be implemented if there ex- 
ist, conditional to the state of control qubit, two different 
pairs of cyclic states of the target qubit. In terms of the 
computational basis {|00), |01), |10), |11)}, where the first 
(second) bit represent the state of the control (target) 
qubit. The unitary operator describing the conditional 
two-qubit gate is given by = diag(U( 7 a iX o\, U( 7 i iX i\), 
under the condition that the control qubit is off reso- 
nance in the manipulation of the target qubit. Here 7* 
(x S ) represents the total phase (the cyclic initial state) 
of the target qubit when the control qubit is in the state 
5(= 0,1). This scheme can be implemented in several 
realistic physical svstems |l2Lll3| . 

Control parameter fluctuation- A simple approach for 
implementing and C/( 2 ' is to use an effective rotating 
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magnetic field to manipulate the state of qubits. In this 
case, the Hamiltonian for single qubit reads 

H = {ujqU x cosujt + ujQtjy smut + lui<j z )/2, (1) 

where u>i — —gfiBi/h (i = 0,1) with g (/i) being the 
gyromagnetic ratio (Bohr magneton), and Bi (i = 1,2) 
acts as an external controllable parameter, and its mag- 
nitude can be experimentally changed. We also use this 
Hamiltonian to manipulate the target qubit in the imple- 
mentation of two qubit gate when the control qubit is off 
resonance. The Hamiltonian of the two-qubit system is 
given by the Hamiltonian Q plus the coupling Hamilto- 
nian acting on the two qubits. In the non-ideal case the 
control fields contain randomly fluctuating components, 
here we assume that uji is flatly distributed in the inter- 
val [(1 — 8i)uji 1 (1 + 8i)oji] with Si a constant, and then 
we numerically calculate the average fidelity of quantum 
gates [A 1 ) and U^ 2 \ But we will assume that u> is not 
affected by random fluctuations; this seems reasonable 
since frequency may be well controlled in a realistic ex- 
periment. 

The average fidelity of a quantum gate we study is 
defined by 

2=1 2=1 

(2) 

where = [cos ^e~ 1 ^/ 2 , sin ^ e ^/ 2 ] T or \ipf) = 

[-sin^-e- l ^/ 2 ,cos^e iw/2 ] T (T denotes the matrix 
transposition) is an input state. 6j G [0, n] and (fj £ 
[0, 2tt] are randomly chosen in our numerical calculation. 
Uid is the ideal quantum gate without any control pa- 
rameter fluctuation and \ip° ut ) is the output state after 
a noisy gate operation when the input state is \ipj n ). In 
the numerical calculation, we randomly choose one input 
state IV'j"), and then calculate the average fidelity of Fj, 
up to satisfactory convergence, for M configurations of 
fluctuations of magnetic fields. After that, we randomly 
choose the next input state and repeat the above calcula- 
tion until deriving the fidelity of this specific input state 
with satisfactory convergence. We repeat N times to cal- 
culate the average fidelity by randomly choosing {ipj 1 ). In 
our numerical calculations below, we get small statistical 
errors when M and N are about several hundreds to one 
thousand. 

Fidelity of single-qubit gates- The Schrodinger equa- 
tion with Hamiltonian can be analytically solved, 
and single qubit gates can be achieved, where 

X = arctan[u;o/(wi — uj)] is the angle between the ini- 
tial state and the symmetric axis of the rotating field. 
The corresponding phases for one cycle are given by jd = 
— tt[u>q+uj-l(u)i— uj)]/uj£1, 7 9 = — 7r[l— (uji— u>)/fl], and 7 = 
— 7r(l + f2/u/) with £1 = ya)g + (uj\ — uj) 2 . We can choose 
any two processes with different values Wg,u^} (j — 
1, 2) satisfying the constrain sin 7! sin 72 sin(x 2 — Xi) 7^ 
to enact two noncommuting single-qubit gates|13j- 



What is remarkable here is that quantum gate im- 
plemented in this way can be varied continuously from 
a standard dynamic gate to a pure geometric gate, by 
changing the external parameters {uj, loq, Wi}. Hence, 
this scheme looks like an ideal model to compare the dif- 
ference between geometric quantum gates and standard 
dynamic gates. It is straightforward to verify from the 
expression of jd that the dynamic phase is zero under the 
condition u> — (uj 2 +u) 2 )/oji. Thus we can obtain purely 
geometric gates by choosing these specific parameters, 
such that 7d = in the whole process. It has been shown 
that x hi can be controlled independently by the 
symmetric axis of the rotating field, and a specific quan- 
tum gate can be realized by a fixed phase 7. For example, 
the Hadamard gate is obtained by 7 — (ni + l/2)-7r and 
X = i n 2 + l/4)7r (with m,2 integers). Therefore, we as- 
sume the total phase is fixed, for concreteness, we choose 
7 = — Ptt with (3 a constant. It is str aightforward to check 
that uj = (w\ ± — 77(^0 + ^1)) / V with rj = 2(3 — (3 2 
guarantee 7 = — f3ir. The obvious requirement of uj re- 
ality implies rjuj 2 < (1 — rf)uj\. On the other hand, the 
dynamic phase is zero under the condition 

uji = uj a \Jrj/(l - 77). (3) 
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FIG. 1: The fidelity and phase in single qubit gates, (a) 
Fidelity, (b) dynamic phase. 
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FIG. 2: The fidelities and phases in single qubit gate for <5o = 
0.1. The values of Si are indicated, (a) Fidelities, (b) dynamic 
and geometric phases. 

The fidelity F(ujq, A) and dynamic phase of this gate 
for typical parameters are shown in Fig.l. Here A is de- 
fined by uj\ = \/r)/(l — r/)ujQ + A. /? — 3/2 is chosen just 
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as an example; we have checked that the main features 
described here are independent on this parameter. We 
plot F(u>o,A) instead of F{u>a,uj{) to guarantee that u> 
(= 2(2wi — y/ojf ~ 3cj 2 )/3) is a real number in the whole 
region. We can see several remarkable features from 
Fig^ where 5q = 8\ = 0.1: (i) the maximum of fidelity 
is along the line described by A = 0, where the dynamic 
phase is zero; (ii) Fidelity is monotonically decreasing by 
increasing A. The change of fidelity is slow at small A, 
but very sharp near A ~ loq, and it is also slow for large 
values of A. Remarkably, the changes of the absolute 
value of dynamic phases with A are just as the same as 
the changes of the average fidelity of gates. Therefore, 
it clearly shows the close relation between the fidelity of 
quantum gates and the component of the dynamic (ge- 
ometric) phase. This feature will be demonstrated also 
in the two-qubit case addressed later. We also checked 
the fidelity for different So and 8\ . The fidelity and quan- 
tum phases as a function of A for ujq = 10 5 are shown 
in Fig [3 We observe that the two main features dis- 
cussed above appear, for the case 8q = 0.1, when 8\ is 
greater than 0.04. However, when Si is less than 0.04, it 
is worth pointing out that the fidelity in the points with 
A = is a local maximum, since there is a dip clearly 
shown nearby A ~ ujq; the largest fidelity appears when 
dynamic phase is dominant. We also numerically com- 
puted the fidelity for fixed but varied Sq (not shown), 
the main features for different Sq are totally similar to 
those in Fig^ The fact that the fidelity is larger when 
A (u>i) is large but with small fluctuation of lo\ can be 
qualitatively explained from the dynamic phase 7<f. The 
deviation of "fd from the noiseless case is dominated by 
the fluctuations of u>i when A is much larger than ojq. 
Therefore the infidelity should be small if the fluctua- 
tions of u>i are very small. 

Fidelity of two-qubit gates- We now numerically com- 
pute the fidelity of a two-qubit gate. We assume that the 
qubit-qubit interaction is given by Hi = Jai 1 ^ cri 2 ^ /2. 
This kind of coupling between qubits can be naturally 
realized, e.g., in quantum computer models with NMR 
and superconducting charge qubits coupled through ca- 
pacitors. When the target qubit is manipulated by a 
rotating field described by Eq.QJ and control qubit is 
off resonance, it is shown that a two-qubit gate U 1 - 2 ^ 



with x = arctan[w /(wi — oj)] and 7 



-tt(1 + Qf/uj) 



can be implemented. Here uif — w\ + (26 — 1) J and 
Q s = \/wq + (wf — lu) 2 , li, uiq, and u>\ are parameters 
for the target qubit. The corresponding phases for one 
cycle are given by 7^ = — 7r[ti> 2 + wf (uif — w)] / ujCI 5 , and 
7* = — 7r[l — (uii — w)/Q s ]. Besides, it is easy to check 
from 7^ = that the geometric two-qubit gates are real- 



ized whenever ui = 2u>\, and uif 



J- 



There is a lot of freedom in choosing parameters to im- 
plement a geometric quantum gate in the present scheme. 
One possible choice is given by u> = u>i + VT+ a 2 Wo 
and J = aujQ, thus the speed of the purely geometric 
gate is of the same order of that of the dynamic gate. 



To see the relation between the fidelity and the quan- 
tum phase, we plot the fidelity and dynamic phase of 
gate U^ 2 ' just when 8 — as a function of the external 
parameters coq and oj\ in Fig. 3. In this case, the dy- 
namic phases change sharply nearby the line described 
by log 10 o;i = log 10 (vl + a 2 wo), where dynamic phases 
are zero, since it is straightforward to find that under the 
condition 



= vT + a 2 u , 



(4) 



7^ is zero. We can see from Fig 01 that the maximum 
of fidelity is along the line where the dynamic phase 
is zero. Moreover, compared with the single-qubit case 
shown in Fig^ the fidelity of the gate shown here de- 
creases quickly, since the dynamic phase changes sharply 
when the parameters do not satisfy Eq.(£j}. Therefore, 
it is clearly shown the close relation between the change 
of fidelity and the change of dynamic component of the 
phase. 




FIG. 3: The fidelity and phase in two-qubit gates for 5q = 
Si = 0.1, 6 = and a = v3- (a) Fidelity, (b) dynamic phase. 
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FIG. 4: The fidelities of two-qubit gates. The curves for dif- 
ferent a are vertically shifted for clarity. Points with zero 
dynamic phases shown in Eq.Jll are denoted by arrows. 

To show in a more clear fashion that the maximum 
of fidelity is along the line described by Eq.QJ and that 
this feature is independent on the parameter a, we nu- 
merically computed the average fidelity of when 
the state of the control qubit is unfixed. The fidelity 
as a function of uq for u\ = 60, 80 = 8± = 0.05 and 
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a = a/3, VS, Vl5, %/35, V143 are plotted in FigH It is 
easy to derive from Eq. that the dynamic phase is zero 
at ujq = 30, 20, 15, 10, 5, and these points are denoted by 
arrows in FigQ] We observe that the maxima of fidelity 
are indeed at the points described by Eq.(@J, where the 
dynamic phases are zero; this property is independent on 
a. 

Comparison with previous results- We would like to 
compare the results here with the previous results in lit- 
erature. It has been shown that the effects of fluctuations 
of the control parameters H7I or decoherence described by 
Lindblad form|lq on nonadiabatic gates are more severe 
than for the standard dynamic gates. We note that in 
Refs.0,0] the geometric gate is implemented by three 
rotations, but only one operation is used to realize a stan- 
dard dynamic gate. Thus a direct comparison is some- 
what not appropriate, and one can not rule out other pos- 
sibilities. In the present model, however, the operations 
for both dynamic gate and geometric gate are totally the 
same, except that the controllable parameters vary con- 
tinuously. In this sense this model looks definitely more 
suitable to assess the difference between geometric gates 
and dynamic gates as far as noise resilience is concerned. 

Before concluding, we would like to point out that the 
scheme studied here is based on nonadiabatic operations. 
Note that to is of the same order of the magnitude as too 
or loi in both one and two-qubit gate operations. This 
implies that the speed of geometric quantum gate here in- 



vestigated is comparable with that of the dynamic quan- 
tum gate. In contrast, the speed of quantum gate based 
on adiabatic Berry's phase is much lower than that of 
gate using dynamic phase, since the adiabatic condition 
requires that both loq and oj\ should be much larger than 
u). Thus the speed constraint required by adiabatic geo- 
metric gates is removed in the present scheme, and the 
results discussed above show the robustness of this nona- 
diabatic geometric gates against noise in the control pa- 
rameters. This considerations suggest that this approach 
to quantum computation should be significant in quan- 
tum information processing. 

Conclusions- In this paper we have numerically com- 
puted the average fidelity of a family of quantum gates, 
which vary continuously from the standard dynamic 
gates to pure geometric gates with the change of exper- 
imentally controllable parameters. We found that the 
maximum of gate fidelity is achieved along the line where 
the dynamical phase is zero. We believe that the results 
presented in this paper provide a first convincing evi- 
dence of the robustness of geometric computation. Our 
predictions can be, in principle, experimentally tested in 
already existing quantum computer prototypes, such as 
NMR, trapped ions and superconducting qubits, etc.. 

This work was supported by the European Union 
project TOPQIP (contract IST-2001-39215). S.L.Z. was 
supported by the NSFC under Grant No. 10204008 and 
the NSF of Guangdong under Grant No.021088. 



[1] See e.g., D.P. DiVincenzo and C. Bennett, Nature 404, 
247 (2000) 

[2] P. W. Shor, Phys. Rev. A 52, R2493 (1995); A. M. 

Steane, Phys. Rev. Lett. 77, 793 (1996). 
[3] L. M. Duan and G. C. Guo, Phys. Rev. Lett. 79, 1953 

(1997); P. Zanardi and M. Rasetti, ibid. 79, 3306 (1997). 
[4] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998); P. 

Zanardi, Phys. Lett. A 258, 77 (1999). 
[5] M. V. Berry, Proc. R. Soc. London A392, 45 (1984). 
[6] Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1953 

(1987). 

[7] S. L. Zhu, Z. D. Wang, and Y. D. Zhang, Phys. Rev. B 

61, 114 (2000). 
[8] P.Zanardi and M.Rasetti, Phys. Lett. A 264, 94 (1999). 
[9] J.A.Jones, V.Vedral, A.Ekert, and G.Castagnoli, Nature 

(London) 403, 869 (2000). 
[10] L.M.Duan, J.I.Cirac, and P.Zoller, Science 292, 1695 

(2001). 

[11] X. B. Wang and M. Keiji, Phys. Rev. Lett. 87, 097901 

(2001) . 

[12] S.L.Zhu and Z.D.Wang, Phys. Rev. Lett. 89, 097902 

(2002) ; Phys. Rev. A 66, 042322 (2002). 

[13] S. L. Zhu and Z. D. Wang, Phys. Rev. A 67, 022319 

(2003) . 



[14] D.Leibfried, B.DeMarco, V.Meyer, D.Lucas, M.Barrett, 
J.Britton, W.M.Itano, B.Jelenkovic, C.Langer, 
TRosenband, and D.J.Wineland, Nature (London) 
422, 412 (2003); S. L. Zhu and Z. D. Wang, Phys. Rev. 
Lett. 91, 187902 (2003). 

[15] P. Solinas, P. Zanardi, and N. Zanghi, quant-ph/0312109 
(2003). 

[16] A. Nazir, T. P. Spiller, and W. J. Munro, Phys. Rev. A 

65, 042303 (2002). 
[17] A. Blais and A. -M. S. Tremblay, Phys. Rev. A 67, 012308 

(2003). 

[18] A. Carollo, I. Fuentes-Guridi, M. F. Santos, and V. Ve- 
dral, Phys. Rev. Lett. 92, 020402 (2004); ibid. 90, 160402 
(2003). 

[19] G. D. Chiara and G. M. Palma, Phys. Rev. Lett. 91, 
090404 (2003). 

[20] Y.Makhlin, G.Schon, and A.Shnirman, Rev. Mod. Phys. 
73, 357 (2001). 

[21] L. M. K. V andersypen and I. L. Chuang, 

|quant-ph /0404064 
[22] J. .1. Garcia-Ripoll and J. I. Cirac, Phys. Rev. Lett. 90. 

127902 (2003). 



